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Charm Dalitz Analyses 
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A review of recent experimental results of Dalitz analyses of charmed meson decays into three-body final states 
is presented. These analyses can help in understanding the strong interaction dynamics leading to the observed 
light mesons spectrum (low mass scalar a, /o{980), ao{980)). A model for the decay amplitude into such 
states is very important for the extraction of the angle 7 of the CKM unitarity triangle. Implications for such 
measurement are discussed. 
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1. Introduction 



A D meson is as a unique "laboratory" to study 
light quark spectroscopy. It has a well defined spin- 
parity — 0~ , constraining the angular momentum 
of the decay products in multibody final states which 
can be analyzed with the Dalitz plot technique [1]. 

Investigations of the low mass scalar mesons can 
be pursued in three-body decays of pseudoscalar D 
mesons giving their large coupling to such states. The 
nature of such low mass scalar states is still under dis- 
cussion _2j , since scalar mesons are difficult to resolve 
experimentally because of their large decay width. 
There are claims for the existence of broad states close 
to threshold such as k(800) and cr(500) On the 
theory side the scalar meson candidates are too nu- 
merous to fit in a single J^^ — qq nonet and 
therefore alternative interpretations are proposed. For 
instance, ao(980) or /o(980) may be 4-quark states due 
to their proximity to the KK threshold Q. 

These hypotheses can be tested through an accurate 
measurement of branching fractions and couplings to 
different final states. In addition, comparison between 
the production of these states in decays of differently 
flavored charmed mesons D°(cu), I?+(cd) and Df{cs) 
can yield new information on their possible quark 
composition. Another benefit of studying charm de- 
cays is that, in some cases, partial wave analyses are 
able to isolate the scalar contribution almost back- 
ground free. 

Results of D° Dalitz analyses can be an in- 
put for extracting the CP-violating phase 7 = 
a.'Tg {-V^^V*f^/V^jV*fj) of the quark mixing matrix by 
exploiting interference structure in the Dalitz plot 
from the decay D°K^ Modeling of the 

Kit and tttt S-wave in D decays is therefore an im- 
portant element in such measurement, since the sys- 
tematic uncertainty on 7 due to the Dalitz model is 
dominated by such components f?]. Model indepen- 
dent approaches using special Dalitz charm analyses 
are discussed and a projected systematic error on 7 in 
future experiments evaluated. 



2. Dalitz analysis formalism 

The amplitudes describing D meson weak-decays 
into three-body final states are dominated by inter- 
mediate resonances that lead to highly non-uniform 
intensity distributions in the available phase space. 

Neglecting CP violation in D meson decays, we de- 
fine the D (D) decay amplitude A (A) in a -D ^ ABC 
Dalitz plot, as: 



A[D ^ ABC] = fDo{mlc,m\c). 
A[D^ BAC] = fooim 



AC' '' 



'■BCl 



(1) 

(2) 



The complex quantum mechanical amplitude / is 
a coherent sum of all relevant quasi-two-body D — s- 
(r — > AB)C resonances ("isobar model" [8]), / = 



''AS' 



and Aj. is the res- 



Y,^are''''-Ar{s). Here s 

onance amplitude. The coefficients and (pr are usu- 
ally obtained from a likelihood fit. The probability 
density function for the signal events is |/|^. Sub- 
modes branching fractions (" fit fractions" ) are defined 
as 



'^j.r '^3^1 1 AjAldm\(jdm\(j 



The fractions fr do not necessarily add up to 1 because 
of interference effects among the amplitudes. 

For well established resonances of the spin-1 {P- 
wave) and spin-2 states, the Breit-Wigner amplitude 
is used 



Abw{s) = Ml{s,p) 



1 



T{s) = F, 



-) 



s ~ iMoF(s) ' 



(3) 
(4) 



where Mq (Fq) is the resonance mass (width) 1], 
L is the angular momentum quantum number, p is 
the momentum of either daughter in the resonance 
rest frame, and po is the value of p when s = Mq. 
The function !Fl is the Blatt-Weisskopf barrier 
factor fl^: J^p = 1, J^i = ^/^/l + Rp^, and T2 

— 1/ \/9 + 3Rp^ + Rp"^, where we take the meson 
radial parameter R is usually set to 1.5 GeV~^ [ll|. 
The spin part of the amplitude, Mli is defined 
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as: 



Ma = Ml = -2 PA-PC, and AI2 = | 



Hpa-PcY 



\PAf-\Pcf] 



M 



D ' 



where Aljy is the 



nominal D mass, and pi is the 3-momentum of 
particle i in the resonance rest frame. 

The Ar{s) parameterization of the scalar /o(980) 
resonance, whose mass, to/q, is close to the KK pro- 
duction threshold, uses the Flatte 03 formula 



^/o(980)("^) = 



Jo 



(5) 



where g/oTr-n- and gj^xTi ^'^^ /o(980) coupling con- 
stants of the resonance to the tttt and KK final states, 
and Pabi-s) = '2pa/\/is) is a phase space factor, cal- 
culated for the decay products momentum, pa, in the 
resonance rest frame. A similar formula is used for 
the ao(980) scalar resonance. 

Different models for the low mass tttt S wave, (called 
a or /o(600) ) are used. In Q a simple spin-0 Breit- 
Wigner is tried. Alternatively a complex pole ampli- 
tude proposed in Ref. [13] can be used 



1 



(6) 



where mc — (0.47 — iO.22) GeV is a pole position 
in the complex s plane estimated from the results of 
several experiments. 

More comprehensive parameterizations of the low 
mass TTTT S wave has been proposed and tested 
[13] [1^ . A K-matrix approach [l6l . \it\ , which gives a 
description of S wave tttt resonances treating the a and 
/o(980) contributions in a unified way has been used 
giving comparable results to the isobar technique [Tsj . 



3. TTTT S-wave. 



3.1. CLEO-c D- 



TX Tl TT 



[19]. 



A study of charged D decay to three charged pions 
has been carried out with the CLEO detector 
This mode has been studied previously by E687 
E691 [21], E791 [1], and FOCUS [18]. 

E791 uses the isobar technique, where each reso- 
nant contribution to the Dalitz plot is modeled as a 
Breit-Wigner amplitude with a complex phase. This 
works well for narrow, well separated resonances, but 
when the resonances are wide and start to overlap, so- 
lutions become ambiguous, and unitarity is violated. 
In contrast, FOCUS uses the K-matrix approach. The 
two techniques give a good description of the observed 
Dalitz plots and agree about the overall contributions 
of the resonances. Both experiments see that about 
half of the fit fraction for this decay is explained by a 
low 7T~^TT~ mass S wave. 



The CLEO analysis utilizes 281 pb~^ of data col- 
lected on the ■0(3770) resonance at y/s ~3773 MeV at 
the Cornell Electron Storage Ring, corresponding to 
a production of about 0.78 x 10^ D^D~ pairs. 
mesons are produced close to the threshold, and are 
thus almost at rest. Events from the decay 
KgTT'^, which has a large rate and contributes to the 
same final state, are isolated with the tt+tt" invariant 
mass even without clearly detached vertexes as in the 
fixed target experiments. 

An isobar model is used to parametrize the signal 
decay where the description of the a from Ref. [11] and 
the Flatte parameterization for the threshold effects 
on the /o(980) Jjj are included. Alternative models 
are also tried and give comparably good fit results 

The TT 7r+7r+ decay tracks are selected with 

requirements on their impact parameters with respect 
to the beam spot. This removes ^--^60% of events with 
Kg — > 7r+7r~ decays. The remaining events from 
Z?+ KgTT^ represent about one third of those se- 
lected for the Dalitz plot. Selection of events from the 
— > 7r^7r+7r+ decay is done with two signal vari- 
ables: AE = Ed~ Ehcmn and mec = V^'Lm-n^ p'd^ 
where -Ebcam is a beam energy, and E^ and po are the 
energy and momentum of the reconstructed D meson 
candidate, respectively. This gives 6991 events in the 
signal box, 2159±18 of these estimated to be back- 
ground. 

The presence of two 7r+ mesons impose a Bose- 
symmetry of the 7r"7r+7r+ final state. The Bose- 
symmetry when interchanging the two same sign 
charged pions is explicitly accounted for in the am- 
plitude parameterization. Dalitz plot is analyzed by 



choosing x = m^iiT'^'K 



J Low 



and y 



m 



')High 



as the independent (x,y) variables. The third variable 
z = m^(7r+7r+) is dependent on x and y through the 
energy-momentum balance equation. 

In the Dalitz plot analysis events in the band 0.2 < 
m^(7r+7r")Low < 0.3 (GeV/c^)^ are excluded which 
is approximately ten times our K'g — > tt^tt" mass 
resolution. This leaves for the Dalitz plot analysis 
4086 events ^2600 of which are signal events. 

CLEOc was able to reproduce the fit results 
E791 [1]. The amplitude normalization and sign con- 
ventions are different from E791, in particular the in- 
clusion of a (TTT contribution gives a fit probability of 
~ 20%. Possible contributions form all known tt+tt^ 
resonances listed in Ref. [2^ were tried, including high 
mass resonances giving asymptotic "tails" at the edge 
of the kinematically allowed region. 

For the /o(980) the Flatte formula, Eq. [5l is used 
with parameters taken from the recent BES II mea- 
surement [llj. For the a a complex pole amplitude, 
Eq. [6l was eventually tried rather than the spin-0 
Brcit-Wigner. 

Table |T] shows the list of surviving contributions 
with their fitted amplitudes and phases, and calcu- 
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Table I Results of the isobar model analysis of the — > 
7r~7r+7r"'" Dalitz plot. For each contribution the relative 
amplitude, phase, and fit fraction is given. The errors are 
statistical and systematic, respectively. 



Mode 


Amplitude (a.u.) 


Phase (°) 


Fit fraction (%) 


p(770)7r+ 


1 (fixed) 


O(fixed) 


20.0±2.3±0.9 


/o (980)71+ 


1.4±0.2±0.2 


12±10±5 


4.1±0.9±0.3 


/2(1270)^+ 


2.1±0.2±0.1 


-123±6±3 


18.2±2.6±0.7 


/o(1370)7r+ 


1.3±0.4±0.2 


-21±15±14 


2.6±1.8±0.6 


/o(1500)7r+ 


1.1±0.3±0.2 


-44±13±16 


3.4±1.0±0.8 


a pole 


3.7±0.3±0.2 


-3±4±2 


41.8±1.4±2.5 



lated fit fractions after a procedure of addition and 
removal of resonances to improve the consistency be- 
tween the model and data. The sum of all fit fractions 
is 90.1%, and the fit probabiUty is ~28% for 90 de- 
grees of freedom. The two projections of the DaUtz 
plot and selected fit components are shown in Fig. [1] 

For contributions that are not significant upper lim- 
its at the 95% confidence level are set. 

The systematic uncertainties, shown in Table [H are 
estimated from numerous fit variations, by adding or 
removing degrees of freedom, changing the event se- 
lection, and varying the efficiency and background pa- 
rameterizations. 

For the poorly established resonances as the a pole, 
their parameters are allowed to fioat and the varia- 
tions of the other fit parameters contribute to the sys- 
tematic errors. The fitted values for the a pole are 
Re{ma) (MeV/c2) = 466±18 and Jm(m<,) (MeV/c^) 
= -223±28. 



X+Y projections I 



Z-projeclion 





selected in a sample of events having at least three 
reconstructed charged tracks with two well identi- 
fied kaons and one pion. The three tracks are fit 
to a common vertex with the constraint they come 
from the beamspot. The decay chain _D*(2112)+ 
D+7 helps in discriminating signal from combinato- 
rial background. Additional requirements based on 
kinematic and geometric information are combined to 
further suppress the background. The final sample 
contains 100850 events with a purity of 95%. An 
unbinned maximum likelihood fit of the Dalitz plot 
(Fig[2]) is performed to extract the relative amplitudes 
and phases of the intermediate resonances as shown 
in Tab.ini The decay is dominated by the (/)(1020)7r+ 
and /o(980)7r+. The /o(980) is parametrized with a 
coupled channel Breit-Wigner [23[ and its contribu- 
tion is large but it is subject to a large systematic 
error due to the poor knowledge of its parameters and 
possible ao(980) contributions that are difficult to dis- 
entangle in the KK projection. Analysis of the angu- 
lar moment distribution confirms such picture with a 
big S-wave-P-wave interference in the KK channel in 
the region of the (/)(1020). On the other hand very 
small activity is present in the K*{892) region sug- 
gesting a small Kir S-wave, and therefore no evidence 
of a k(800). 





m-(K*K")(GcV7c'b 



m-(K* K") (GeVV) 




05 1 1.5 2 

m-(K II*) (GeV'^/c*) 




m'CKV) (GcV V) 



Figure 1: Projection of the Dalitz plot onto the 
m^(7r"'"7r~) axis (two combinations per D+ candidate) for 
CLEO-c data (points) and isobar model fit (histograms) 
showing the various components (left) .Projection of the 
Dalitz plot onto the m^(7r+7r+) axis (right). 



Figure 2: Dalitz plot of Ds ~* tv+K+R- 



3.3. BaBar 



3.2. BaBar D, 



K+K^n' 



BaBar analyzed 240 fb ^ taken at the center of 
mass energies near the T{4S) resonance. Events are 



The data sample used in the BaBar 
K^K~K~^ analysis consists of 91.5 fb~^ recorded with 
the BaBar detector at the SLAG PEP-II storage rings 
24| . The PEP-II facility operates nominally at the 



fpcp07_252 



4 



Flavor Physics and CP Violation Conference, Bled, 2007 



Table II Results of the isobar model analysis of the Ds — ^ K K n Dalitz plot. For each contribution the relative 
amplitude, phase, and fit fraction is given. The errors are statistical and systematic, respectively. 



Mode 


Amplitude (a.u.) 


Phase (°) 


Fit fraction (%) 


K*{m2)K+ 


1 (fixed) 


O(fixed) 


48.7 ±0.2 ±1.6 


<?!>(1020)7r+ 


1.081 ± 0.006 ± 0.049 


2.56 ±0.02 ±0.38 


37.9 ±0.2 ±1.8 


/o(980)7r+ 


4.6 ±0.1 ±1.6 


-1.04 ±0.04 ±0.48 


35 ± 1 ± 14 




1.07 ±0.06 ±0.73 


-1.37 ±0.05 ±0.81 


2.0 ±0.2 ±3.3 


/o(1710)^+ 


0.83 ±0.02 ±0.18 


-2.11 ±0.05 ±0.42 


2.0 ±0.1 ±1.0 


/o(1370)7r+ 


1.74 ± 0.09 ± 1.05 


-2.6 ±0.1 ±1.1 


6.3 ±0.6 ±4.8 




0.43 ± 0.05 ± 0.34 


-2.5 ±0.1 ±0.3 


0.17 ±0.05 ±0.30 


/2(1270)7r+ 


0.40 ± 0.04 ± 0.35 


0.3 ±0.2 ±0.5 


0.18 ±0.03 ±0.40 


Sum 


132 ± 1 ± 16 
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Figure 3: Dalitz plot of D° ^ Ti'^K+K-. 



ao(980)+ can be observed on the lower right edge of 
the Dalitz plot. This contribution, however, does not 
overlap with the 0(1020) region and this allows the 
K^K~ scalar and vector components to be separated 
using a partial wave analysis in the low mass 
region. 

The helicity angle, 6*^, is then defined as the angle 
between the K+ for D° (or K' for W) in the K+K' 
rest frame and the K'^K~ direction in the (or 
K^) rest frame. The K^K^ mass distribution has 
been modified by weighting each candidate by the 
spherical harmonic {cos 9k) (L=0-4) divided by its 
(Dalitz-plot-dependent) fitted efficiency, ft is found 
that all the ) moments are small or consistent with 
zero, except for (Y^), (Y^) and (Y^). 

In order to interpret these distributions a simple 
partial wave analysis has been performed, involving 
only S- and P-wave amplitudes. This results in the 
following set of equations (25j : 



T(45) resonance, providing collisions of 9.0 GeV elec- 
trons on 3.1 GeV positrons. The data set includes 
82 fb~^ collected in this configuration (on-resonance) 
and 9.6 fb~^ collected at a cm. energy 40 MeV below 
the T{AS) resonance (off-resonance). 

Selecting events within ±2cr of the fitted D'^ mass 
value, a signal fraction of 97.3% is obtained for the 
12540 events selected. The Dalitz plot for these 
1?° K^K^K" candidates is shown in Fig. [3l In 
the K~^K~ threshold region, a strong (/)(1020) signal 
is observed, together with a rather broad structure. 
A large asymmetry with respect to the K'^K"^ axis 
can also be seen in the vicinity of the 0(1020) signal, 
which is most probably the result of interference be- 
tween S and P-wave amplitude contributions to the 
K+K- system. The /o(980) and ao(980) S'-wave res- 
onances are, in fact, just below the K^K^ threshold, 
and might be expected to contribute in the vicinity of 
(/)(1020). An accumulation of events due to a charged 



V^{Y^) ^S^ + P^ 
VI^{Y,°)^2\S\\P\coscl)sP (3) 

where S and P are proportional to the size of the 5- 
and P-wave contributions and (psp is their relative 
phase. Under these assumptions, the (^2") moment 
is proportional to P^ so that it is natural that the 
0(1020) appears free of background, as is observed. 
This distribution has been fit using the following rel- 
ativistic P-wave Breit-Wigner, yielding the following 
parameters: 

= 1019.63 ± 0.07, = 4.28 ± 0.13 MeV/c^ 

in agreement with PDG values (statistical errors 
only). 
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The above system of equations can be solved di- 
rectly for S^, P^ and cos 4>sp and corrected for phase 
space distribution. The phase space corrected spectra 
are shown in Fig. 01 

The distributions have been fitted using a model 
with (/)(1020) for the P-wave , a scalar contribution 
in the K'^K^ mass projection entirely due to the 
ao(980)*', K'^K^ mass distribution is entirely due to 
ao(980)+ and the cos ^sp described with BW models. 

The ao(980) scalar resonance has a mass very close 
to the KK threshold and decays mostly to rjir. It has 
been described by a coupled channel Breit Wigner of 
the form: 



BWch{ao){m) = 



9kk 



(5) 



where pirn) — 2q/m while grjTv and g^^^ describe the 
ao(980) couphngs to the r?7r and KK systems respec- 
tively. 

The best measurements of the ao(980) parameters 
come from the Crystal Barrel experiment j2|||, in 
pp annihilations, with a value of giix — 329 ± 27 
(MeV)^/^. mo and have been fixed to the Crystal 
Barrel measurements, but g^K^ the other hand, 
has been fit (stat only) g^x ^ 464 ± 29 (MeV)^/^ 

The determination of gj^}^ has been redone in a 
complete Dalitz plot analysis with an evaluation of 
the systematic error. The fit produces a reasonable 
representation of the data for all of the projections. 
The computed on the Dalitz plot gives a value 
of X /NDF=%%ilnA. The sum of the fractions is 
130.7 ± 2.2 ± 8.4%. The regions of higher are 
distributed rather uniformly on the Dalitz plot. At- 
tempts to improve the fit quality by including other 
scalar amplitudes caused the fit to diverge, producing 
a sum of fractions well above 200% along with small 
improvements of the fit quality. 

The final fit results showing fractions, ampli- 
tudes and phases are summarised in Table IIIII For 
7?o/o(980) and ii'+ao(980)- (DCS), being consistent 
with zero, only the fractions have been tabulated. For 
the Dalitz plot analysis the /o(980) contribution is 
found to be consistent with zero, 

A test has been performed by leaving g^K as a free 
parameter in the Dalitz plot analysis, the resulting 
central value of g^x being 



g^j^ = 473 ± 29 (stat.) ± 40 (syst.)(MeV) 



1/2 



This value differs significantly from the Crystal Barrel 
measurement. An improvement of this measurement 
can be foreseen by adding data from the ao(980) — > rjir 
decay mode such as D° K^rjir'^. 

It must be noticed that reliable estimate of the ex- 
pected contribution of the /o(980) in D° K^K+K' 
decay is not possible until more accurate measure- 
ments of the /o(980) parameters and couplings be- 
come available. This can be performed, for example. 



by using high statistics samples of Df KKtt^ and 
Df — > Tr+vr+vr^ decays. 
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Figure 4: Results from the K^K^ Partial Wave Analysis 
corrected for phase space, (a) P-wave strength, (b) S- 
wave strength, (c) m{K^ K'^) distribution, (d) cos^sp in 
the 0(1020) region, (e) (jjsp in the threshold region after 
having subtracted the fitted 0(1020) phase motion shown 
in (d). The lines correspond to the fit described in the 
text. 



4. Kit S-wave 

The K^n'^ systems from the decay D*' 
K~K^'!r^ can provide information on the Kir S- 
wave (spin-0) amplitude in the mass range 0.6-1.4 
GeV/c^, and hence on the possible existence of the 
k(800), reported to date only in the neutral state 
(k" ^ i4r~7r+) 3]. If the k has isospin 1/2, it should 
be observable also in the charged states. 



4.1. BaBar 



K+K- 



BaBar analyzed 385 fb~^ of e+e^ collision data 
and reconstructed the decays £'°7r+ with 

K~K'^Tr^ [23|. Requirements on the center-of- 
mass momentum of the D'^ candidate and on \m]j,+ — 
m£)o| yields in the signal region, 1855 < mjjo < 
1875 MeV/c^ 11278 ± 110 signal events with a purity 
of about 98.1%. 

For D'^ decays to K^tt'^ S-wa.ve states, three ampli- 
tude models have been considered. One model uses 
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Table III Results from the Dalitz plot analysis of D K K K . The fits have been performed using the value of 
Qjiii — 464 (MeV)'^'''^ resulting from the partial wave analysis. 



Final state 


Amplitude 


Phase (radians) 


Fraction (%) 




1. 


0. 


66.4 ± 1.6 ± 7.0 


ZV(1020) 


0.437 ± 0.006 ± 0.060 


1.91 ± 0.02 ± 0.10 


45.9 ± 0.7 ± 0.7 


X-ao(980)+ 


0.460 ± 0.017 ± 0.056 


3.59 ± 0.05 ± 0.20 


13.4 ± 1.1 ± 3.7 


:^°/o(1400) 


0.435 ± 0.033 ± 0.162 


-2.63 ± 0.10 ± .71 


3.8 ± 0.7 ± 2.3 


:^°/o(980) 






0.4 ± 0.2 ± 0.8 








0.8 ± 0.3 ± 0.8 


Sum 






130.7 ± 2.2 ± 8.4 



Table IV The results obtained from the D K~ K'^n Dalitz plot fit. We define amplitude coefficients, ar and 
(j)r, relative to those of the 7('*(892)"^. The errors are statistical and systematic, respectively. We show the ao(980) 
contribution, when it is included in place of the /o(980), in square brackets. We denote the K-k S-wave states here by 



State 


Amplitude, Ur 


Phase, (l>r (°) 


Fraction, fr (%) 


i^"*(892) + 


1.0 (fixed) 


0.0 (fixed) 


45.2±0.8±0.6 


7^*(1410)+ 


2.29±0.37±0.20 


86.7±12.0±9.6 


3.7±1.1±1.1 


K+n°{S) 


1.76±0.36±0.18 


-179.8±21.3±12.3 


16.3±3.4±2.1 


0(1020) 


0.69±0.01±0.02 


-20.7±13.6±9.3 


19.3±0.6±0.4 


/o(980) 


0.51±0.07±0.04 


-177.5±13.7±8.6 


6.7±1.4±1.2 


[ao(980)°] 


[0.48±0.08±0.04] 


[-154.0±14.1±8.6] 


[6.0±1.8±1.2] 


/2(1525) 


1.11±0.38±0.28 


-18.7±19.3±13.6 


0.08±0.04±0.05 


7t"*(892)- 


0.601±0.011±0.011 


-37.0±1.9±2.2 


16.0±0.8±0.6 


7^*(1410)- 


2.63±0.51±0.47 


-172.0±6.6±6.2 


4.8±1.8±1.2 




0.70±0.27±0.24 


133.2±22.5±25.2 


2.7±1.4±0.8 



the LASS amplitude for if 7r+ - 
tering 28], 

,5(s)=cot-if— + ^Wcot-i 
\pa 2 / 



K 7r+ elastic scat- 



(7) 



,(8) 



where Mq {Tq) refers to the K^iMBO) mass (width), 
a = 1.95 ± 0.09 GeV-^c, and 6 = 1.76 ± 0.36 GeV'^c. 
The unitary nature of Eq. [7] provides a good descrip- 
tion of the amplitude up to 1.45 GeV/c^ (i.e., Kij' 
threshold) . In Eq. [H the first term is a non-resonant 
contribution defined by a scattering length a and an 
effective range b, and the second term represents the 
i4rQ(1430) resonance. The phase space factor \/s/p 
converts the scattering amplitude to the invariant am- 
plitude. 

A second model uses the E791 results for the K^n^ 
S-wsLve amplitude from an energy-independent partial 
wave analysis in the decay D'^ K~Tr~^iT^ (soj . 
The third model uses a coherent sum of a uniform 
non-resonant term, and Breit-Wigner terms for the 
k(800) and i4ro(143G) resonances. 



The Z?°decay to a K^K^ S-wave state is described 
by a coupled-channel Breit-Wigner amplitude for the 
/o(980) and ao(98G) resonances, with their respective 
couplings to TTvr, KK and rjir, KK final states (l^ . 



A 



fa [ao 



Ml. 



Ml 



i{9i PtttvI^tt] + gl Pkk) ' 



(9) 



Several models are used incorporating various com- 
binations of intermediate states. In each fit, the 
if* (892)+ is included and the complex amplitude co- 
efiicients of other states relative to it is measured. 

The LASS Kn S-wave amplitude gives the best 
agreement with data and it is uses it in the nominal 
fits. The Kit S-wave modeled by the combination of 
k(8GG) (with parameters taken from Ref. 0]), a non- 
resonant term and ii'o(1430) has a smaller fit proba- 
bility (x^ probability < 5%). The best fit with this 
model (x^ probability 13%) yields a charged k of mass 
(87G ± 30) MeV/c^, and width (15G ± 20) MeV/c^, 
significantly different from those reported in Ref. [1^ 
for the neutral state. This does not support the hy- 
pothesis that production of a charged, scalar k is be- 
ing observed. The E791 amphtude [301 describes the 
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Figure 5: (Color online) Dalitz plot for D° K-K+tv° 
data (a), and the corresponding squared invariant mass 
projections (b-d). The three-body invariant mass of the 
D'^ candidate is constrained to the nominal value. In 
plots (b-d), the dots (with error bars, black) are data 
points and the solid lines (blue) correspond to the best 
isobar fit models. 



data well, except near threshold (x^ probabihty 23%). 
Analysis of moments of cos 9h confirms little varia- 
tion in S-wave phase up to about 1.02-1.03 GeV/c^ 
and matched the behaviour obtained with the isobar 
model. 

The results of the best fit are summarized in Ta- 
ble IfVI Neglecting CP violation, the strong phase 

difference. So, between the D and I?°decays to 
K*{892)^K^ state and their amplitude ratio, r_D, are 
given by 



5. Dalitz model and CKM 7 extraction 



Assuming no CP asymmetry in D decays the 



decay chain rate 



r^(m^,m^) can be written as 



T^{ml,ml)^\AD^? +rl\AD±? + (H) 

where and rn^ are the squared invariant masses of 
the K'^tt^ and K'^tt^ combinations respectively from 



with Ad~ 



(If 



the decay, and Ad^ 

{Ad+) the amplitude of the D° -> K^n' 
iiTgTr+Tr") decay. In Eq. (11) the following definitions 
are used, = r^ cos{Sb T 7) and j/zp = sin((5s T 
7). Here, is the magnitude of the ratio of the 

amplitudes A{B- ^ ifx*-) and A{B- D°K*-) 
and 5b is their relative strong phase. 




m;„. (GeV7c') 



iSo aDO^K'--K+ ifS 

roe — e ^ 



(10) 



Combining the results of models I and II, we find 5d 
= -35.5° ± 1.9° (stat) ±2.2° (syst) and ro = 0.599 ± 
0.013 (stat) ± O.OII (syst). These results are con- 
sistent with the previous measurements j3l| . 5d = 
-28° ± 8° (stat) ±11° (syst) and ro = 0.52 ± 0.05 
(stat) ± 0.04 (syst). 

The measurement of r^, and is a prerequisite to 
extract the CKM angle 7 from the analysis of B^ 
D^K^ decays [3^], where the symbol indicates ei- 
ther a D^or a D meson decaying into a CP-eigenstate 
as K~K~^Tr° or K'^Tr~Tr~^ as we will see in more detail 
in the next Section. 



Figure 6: (a) The D° KsTv~n+ Dalitz distribution 
from D*~ — » D°Ti-~ events, and projections on (b) = 
'^KOn+' (c) = "1^0 TT-' a^iid (d) m^+^_. D° ^ 

KsTv^n~ f rom D*+ D°Tr+ events are also included. 
The curves are the reference model fit projections. 



Once the decay amplitude Ad is known, the Dalitz 
plot distributions for D'^ from B^ and decays 
can be simultaneously fitted to r_(m^,TO^) and 



r+(TO 



^ , TO^ ) as given by Eq. 



(11), respectively, and 



the angle 7 can be extracted. 

Since the measurement of 7 arises from the interfer- 
ence term in Eq. (11), the uncertainty in the knowl- 
edge of the complex form of can lead to a system- 
atic uncertainty. 
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Table V Complex amplitudes are'*^"" and fit fractions of the different components (A'gTr" and K^tv'^ resonances, and 
7r"'"7r~ poles) obtained from the fit of the K''^tv^tv^ Dalitz distribution from Z)*+ D^n^ events. Errors are 

statistical only. Masses and widths of all resonances are taken from ^25], while the pole masses and scattering data are 
from [sj. The fit fraction is defined for the resonance terms (tttt S-wave term) as the integral of a'^\Ar{rn^-,m\)^ over 
the Dalitz plane divided by the integral of |^i5 (mi , m^)|^. The sum of fit fractions is 1.16. 



Component 


i?e{a,.e^'*'-} 


/mfare"*"" } 


Fit fraction (%) 


if * (892)- 


-1.159 ± 0.022 


1.361 ±0.020 


58.9 


7fo*(1430)- 


2.482 ± 0.075 


-0.653 ±0.073 


9.1 




0.852 ± 0.042 


-0.729 ±0.051 


3.1 


7f*(1410)- 


-0.402 ± 0.076 


0.050 ± 0.072 


0.2 


a:*(1680)- 


-1.00 ±0.29 


1.69 ±0.28 


1.4 


if*(892)+ 


0.133 ±0.008 


-0.132 ±0.007 


0.7 


7^0* (1430)+ 


0.375 ± 0.060 


-0.143 ±0.066 


0.2 


A'2*(1430)+ 


0.088 ± 0.037 


-0.057 ±0.038 


0.0 


p(770) 


1 (fixed) 


(fixed) 


22.3 


w(782) 


-0.0182 ±0.0019 


0.0367 ± 0.0014 


0.6 


/2(1270) 


0.787 ± 0.039 


-0.397 ±0.049 


2.7 


p(1450) 


0.405 ± 0.079 


-0.458 ±0.116 


0.3 


/3i 


-3.78 ±0.13 


1.23 ±0.16 






9.55 ±0.20 


3.43 ± 0.40 




Pa 


12.97 ±0.67 


1.27 ±0.66 




/■prod 
J 11 


-10.22 ± 0.32 


-6.35 ±0.39 




sum of TT+TT^ S-wavc 






16.2 



Two different models describing the D K^tt^tt^ 
decay have been used in this analysis. The first model 
(also referred to as Breit-Wigner model) [111 expresses 
as a sum of two-body decay-matrix elements and a 
non-resonant contribution. In the second model (here- 
after referred to as the tttt S-wave K-matrix model) the 
treatment of the tttt S-wave states in — > ifg7r~7r+ 
uses a K-matrix formalism [H, |T3] to account for the 
non-trivial dynamics due to the presence of broad and 
overlapping resonances. The two models have been 
obtained using a high statistics flavor tagged D'^ sam- 
ple (!?*+ D^Tif) selected from e+e^ — > cc events 
recorded by BaBar. 

In the Breit-Wigner model a set of several two- 
body amplitudes is used, including five Cabibbo- 
allowed amplitudes: i4:*(892)+7r-, X*(1410)+7r-, 
K^{U30)+t:-, ii:2*(1430)+7r- and K*{1680)+t:~ , 
their doubly Cabibbo-suppressed partners, and eight 
channels with a Kg and a tttt resonance: p, u), /o(980), 
/2(1270), /o(1370), /9(1450), ai and cts • The Breit- 
Wigner masses and widths of the scalars cti and CT2 are 
left unconstrained, while the parameters of the other 
resonances are taken to be the same as in [s^]- The 
parameters of the a resonances obtained in the fit are 
as follows: = 519 ± 6 MeV/c^, F^, = 454 ± 12 

MeV/c2, M^^ = 1050 ± 8 MeV/c^ and F,,^ = 101 ± 7 
MeV/c^ (the errors are statistical only). The alter- 
native model is based on a fit to scattering data (K- 



matrix [3J|) used to parametrize the tttt S-wave com- 
ponent. This variation is used to estimate the model 
systematic uncertainty on 7 since it gives an equally 
good fit to data. 

The error due to the resonance model can be 
avoided by using the model-independent 7 measure- 
ment proposed in Q. In this approach, the Dalitz 
plot is partitioned in bins symmetric with respect to 
the 7r+7r~ axis. Counting the number of events in 
such bins from entangled D decay samples, in ad- 
dition to the already utilized fiavour-tagged D de- 
cay samples, can determine the strong phase varia- 
tion over the Dalitz plot. For this the data of a r- 
charm factory is needed. Useful samples consist of 

^^(3770) D°D events where one of the D mesons 
decays into a CP eigenstate (such as K^K~ or i^gw), 
while the D meson going in the opposite direction 
decays into K^gir^ir^ . Using also a similar sample 
where both mesons from the ■(/;(3770) decay into the 
K^TT^'K^ state provides enough information to mea- 
sure all the needed hadronic parameters in D decay up 
to one overall discrete ambiguity (this can be resolved 
using a Breit-Wigner model). CLEO-c showed that 
with the current integrated luminosity of 280 pb~^ 
at the V'(3770) resonance, these samples are already 
available. 

With the luminosity of 750 pb-\ that CLEO-c 
should get at the end of its operation, the samples 
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will be respectively about 1000 and 2000 events. Using 
these two samples with a binned analysis and assum- 
ing rB= 0.1, a 4° precision on <^3 could be obtained 



6. Conclusions. 

Charm meson multi-body decays are crucial to de- 
termine light strong interaction bound states. The 
nature of such mesons is still unclear, but more infor- 
mation is emerging from high statistics Dalitz analy- 
sis of D decays. In the future multi-channles analyses 
may be the way to go to identify underline structure 
of the light mesons. For instance a measurement of 
the couplings of the S-wave in various decays can 
help in interpreting the /o (980) as two di-quark bound 
states [3^ . Determining the decay dynamic of charm 
mesons is relevant for method to extract the CKM 
angle 7 in B decays as B'^ D'^K'^ . The effect of 
the knowledge of the strong phase variation in charm 
meson decay translates into a model systematic error 
on the 7 value. Model dependence can be removed if 
special sample of D meson charm decays in quantum- 
coherent states will be available, bringing down the 
model error on 7 to few degrees. 
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